Interior design of a two-dimensional semiclassical black hole: Quantum transition 

across the singularity 
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We study the internal structure of a two-dimensional dilatonic evaporating black hole, based on 
the CGHS model. At the semiclassical level, a (weak) spacelike singularity was previously found to 
develop inside the black hole. We employ here a simplified quantum formulation of spacetime dy- 
namics in the neighborhood of this singularity, using a minisuperspace-like approach. Quantum evo- 
lution is found to be regular and well-defined at the semiclassical singularity. A well-localized initial 
wave-packet propagating towards the singularity bounces off the latter and retains its well-localized 
form. Our simplified quantum treatment thus suggests that spacetime may extend semiclassically 
beyond the singularity, and also signifies the specific extension. 

PACS numbers: 



I. INTRODUCTION 

The inevitable occurrence of spacetime singularities in- 
side black holes (BHs) constitutes one of the greatest 
challenges in understanding the final fate of gravitational 
collapse. A spacetime singularity is typically a region in 
spacetime where the curvature blows up. The presence of 
spacetime singularities inside BHs has been established 
by several theorems [l| , and is also demonstrated by nu- 
merous examples of exact BH solutions. 

Such curvature singularities are often regarded in the 
literature (at least formally) as the boundary of space- 
time. It is not obvious that this assertion is fully justified 
from the physical view-point (especially in the case of 
weak singularities; sec below). This formal point of view 
is nevertheless pragmatic, because even if spacetime does 
physically extend beyond the singularity, in most cases it 
is unclear what this extension is: The singularity makes 
the mathematical extension ambiguous, as the field equa- 
tions become ill-defined there. One may hope, however, 
that once the full theory of Quantum Gravity is formu- 
lated, it will somehow resolve the spacetime singularities 
and thereby clarify how physics extends beyond them. 

The remarkable revelation of Hawking radiation |2| and 
black-hole evaporation brought new and intense interest 
into black-hole structure, and particularly into the as- 
sociated information puzzle. The latter appears to be 
intimately related to the presence of a singularity in- 
side the BH: As long as the hypersurface of singularity 
is treated as a boundary of spacetime, all bits of infor- 
mation which encounter this hypersurface are quite in- 
evitably "destroyed" there. For this reason one may nat- 
urally hope that the information puzzle will eventually 
be resolved by Quantum Gravity eliminating spacetime 
singularities. 

Two decades ago Callan, Giddings, Harvey and Stro- 
minger (CGHS) [3| introduced a toy model for investigat- 
ing the formation and evaporation of black holes. Their 
toy model consists of two-dimensional gravity coupled to 
a dilaton field. In this two-dimensional framework the 
renormalizcd stress-energy tensor T a p may be expressed 



in a simple explicit form, for any prescribed background 
metric. This property makes the CGHS model a power- 
ful tool for analyzing BH evaporation. At the classical 
level (ignoring the semiclassical corrections), one finds a 
one-parameter family of vacuum solutions, analogous to 
the Schwarzschild solution in four dimensions. Once the 
semiclassical effects arc added, the BH emits a thermal 
Hawking radiation and evaporates. 

Originally it was hoped that the CGHS model of BH 
formation and evaporation will be free of any singularity, 
which would make it an ideal model for addressing the 
information puzzle. However, Russo, Susskind and Thor- 
lacius j4| soon demonstrated that a spacclikc curvature 
singularity actually forms inside the semiclassical CGHS 
BH. This singularity was recently analyzed in some detail 
@ , based on the homogeneous approximation. Quite re- 
markably, the singularity turns out to be weak in Tipler's 
terminology [6[ p}. This means that extended physical 
objects arrive at the singular hypersurface intact (as op- 
posed to a strong singularity where the object is totally 
torn apart). What diverges at the CGHS singularity is 
the rate of change of the metric and the dilaton (e.g. in 
terms of the proper time along a timelike geodesic). In 
other words, the dilaton and metric are continuous but 
non-smooth at the singularity. 

Owing to the weakness of the semiclassical CGHS sin- 
gularity, one may be tempted to assume that semiclas- 
sical spacetime will somehow extend beyond the singu- 
lar hypersurface. This leads one to consider the issue of 
existence and uniqueness of solutions to the semiclassi- 
cal CGHS field equations, in the spacetime region laying 
at the future side of the spacelike singularity — solutions 
which continuously match the data coming from the past 
side. It is not difficult to see that existence is not a prob- 
lem, but uniqueness fails here because we had to give up 
on smoothness. In fact, there is an infinite set of such 
continuous semiclassical extensions beyond the singular- 
ity. Semiclassical theory alone thus fails to predict the 
correct extension — even if such an extension does exist in 
reality One may hope, however, that a fully-quantized 
version of the CGHS model will address this problem 



(e.g. by resolving the singularity). It should tell us if 
spacetime indeed extends semiclassically beyond the sin- 
gularity, and if so — what the correct extension is. 

Motivated primarily by the BH information puzzle, 
Ashtckar, Taveras and Varadarajan (ATV) Q recently 
proposed a quantized version of the CGHS model. In 
their model (roughly speaking) the dilaton and the metric 
conformal factor are elevated to quantum fields defined 
on an -R 2 manifold. These quantum fields are represented 
by operators, within the Heisenberg picture. The field 
equations thus turn into a nonlinear system of operator 
partial differential equations, with initial data specified at 
past null infinity. A key question underlying this model is 
whether the quantum evolution, formulated in this way, 
retains its regularity at the (would-be) semiclassical sin- 
gularity. It is conceivable that regularity is preserved at 
the quantum level, though this still needs to be shown. If 
regularity is indeed preserved, then in principle quantum 
evolution could be followed within this model all the way 
from the collapse to full evaporation. In particular this 
would clarify the state of the quantum fields accessible to 
an observer at future null infinity, thereby resolving the 
information puzzle. Unfortunately the nonlinear opera- 
tor equations underlying the ATV model are extremely 
hard to solve. To overcome this difficulty, ATV intro- 
duced two useful approximates: Bootstrapping, which 
may be applied to the early stage of evaporation; and 
the mean-field approximation, which presumably holds 
near future null infinity. 

In the ATV scenario, the singularity which was present 
in the semiclassical picture is now replaced by a region 
of strong quantum fluctuations. Quantum dynamics is 
presumably regular there, which allows for a well-defined 
evolution of (quantum) spacetime across this region. A- 
priory there are several possibilities: (i) After a narrow 
strongly-fluctuating transition region (marking the locus 
of the semiclassical spacclike singularity), the fluctua- 
tions weaken and spacetime retains its semiclassical char- 
acter; or, alternatively, (ii) quantum fluctuations never 
decline, and spacetime remains strongly-fluctuating in 
the entire causal future of the semiclassical singular- 
ity. Intermediate options are also possible: For example, 
spacetime may become semiclassical again only in the 
neighborhood of future null infinity. 

Our main goal in this paper is to explore the quantum 
transition across the semiclassical spacelike singularity, 
and to reveal the nature of the spacetime region at the 
immediate neighborhood of this singularity. Particularly 
we would like to find out if spacetime indeed becomes 
semiclassical again after crossing the singularity [option 
(i) above], and if so — what is the semiclassical solution 
which takes place there. Unfortunately, none of the ap- 
proximate solutions constructed by ATV apply to this 
portion of the manifold (and the operator field equations 
are extremely difficult to solve). In order to circumvent 
this difficulty, in this paper we shall apply one further 
simplification to the quantum ATV dynamics, and treat 
it at the "minisuperspace" level — that is, we shall ignore 



all degrees of freedom associated with inhomogeneities. 

Our motivation comes from the observation that the 
semiclassical dynamics inside a macroscopic evaporating 
CGHS BH — and in particular the approach to the spacc- 
like singularity — may well be approximated (locally) by 
homogeneous solutions of the CGHS field equations (see 
discussion in Q). In this sense the BH interior resem- 
bles an approximately-homogeneous cosmology, which is 
amendable to a quantum treatment at the minisuper- 
space level. Our quantization strategy thus proceeds as 
follows: (i) We employ the homogenous approximation 
for the BH interior already at the semiclassical level, eras- 
ing any reference to spatial dependence. Correspond- 
ingly we regard our variables (metric and dilaton) as 
mere time-dependent variables, rather than fields; (ii) 
We explore the leading-order asymptotic behavior of our 
variables on approaching the spacelike singularity, (iii) 
We derive an effective Hamiltonian which reproduces this 
leading-order dynamics near the singularity. [Note that 
stages (i-iii), which are implemented at the semiclassical 
level, were already carried out in Ref. Q.] (iv) We now 
quantize our variables, but we do this in a way one quan- 
tizes time-dependent variables (rather than fields). Thus, 
the dynamics is now represented by a time-dependent 
wave- function "J, and the effective Hamiltonian becomes 
a differential operator which determines the time evolu- 
tion of \1/ via a Schrodinger-type equation. 

A key objective of this paper is to explore the prop- 
erties of this simplified quantum system. Wc first ob- 
tain the stationary eigenf unctions, construct from them 
a well- localized initial wave packet, and follow its time 
evolution after it hits the singularity. 

Summary of main results - A key result in our analysis 
is that the quantum evolution is regular and hence well- 
defined. This is not a trivial result, since the Hamiltonian 
operator H itself is in fact singular at the locus of the 
semiclassical singularity. Yet, all its eigenfunctions are 
perfectly regular there, yielding a unique regular quan- 
tum evolution across the semiclassical singularity. 

Our initial conditions correspond to a well-localized 
wave packet which propagates towards the singularity. 
At this initial stage the wave packet follows a semiclas- 
sical orbit, as expected. As the wave packet approaches 
the singularity, it slightly spreads and develops numer- 
ous wiggles. Subsequently, the wave packet bounces 
off the singularity, and retains its original well-localized 
character. In this final stage, too, the wave packet is 
found to follow a semiclassical orbit — a particular mem- 
ber of the infinite set of possible (continuous) semiclas- 
sical solutions beyond the singularity. Thus, this quan- 
tum model determines the specific semiclassical extension 
which takes place to the future of the semiclassical sin- 
gularity. The extension of the dilaton turns out to be (at 
the leading order near the singularity) a time-reflection 
of the evolution prior to the singularity. However, the 
extension of the other degree of freedom (the metric con- 
formal factor) does not admit the same time-reflection 
symmetry. 



Concluding, our quantum model suggests that in an 
evaporating CGHS BH, spacetime physics extends be- 
yond the spacelikc singularity by means of a well-defined 
semiclassical spacetime. It remains to explore the de- 
tailed structure of this new patch of spacetime. In this 
paper we only obtained its leading-order asymptotic be- 
havior just beyond the singularity. The extension beyond 
leading order is discussed in sec. VII. 

We must recall, however, that the quantum model con- 
sidered here is a very simplified one. To obtain more 
firm and trustable results, one must resort to a fuller 
quantum-gravitational theory of spacetime dynamics. 

The paper is organized as follows: In section II we 
briefly review the main results of the semiclassical anal- 
ysis carried out in Ref. [5|, which are relevant to the 
present work. In particular this includes the leading- 
order semiclassical asymptotic behavior near the singu- 
larity, and the corresponding effective Hamiltonian. 

Section III describes our quantization procedure, 
and also presents the eigenfunctions of the stationary 
Schrodingcr equation. Section IV will be devoted to the 
construction of a well-localized initial wave-packet, while 
Sec. V will discuss the time evolution of the wave packet 
and its properties as it hits the singularity and afterward. 
In Sec. VI we present the (leading-order) semiclassical 
extension of spacetime beyond the singularity, as deter- 
mined from the orbit of the well-localized wave packet 
after it bounces off the singularity. Finally, Sec. VII is 
dedicated to discussion of our results and their signifi- 
cance. 

Throughout this paper we use general-relativistic units 
c=G = l. 



II. BACKGROUND 

In Ref. |5| we investigated the interior of a macroscopic 
evaporating CGHS BH, and particularly the singularity 
which develops inside the BH. Here we shall briefly out- 
line the main results, and establish the notation relevant 
to the quantum analysis below. 



For convenience we introduce new variables: R = er 2< ^ 
(which was denoted R in Q) and S = 2(p — 4>). The 
evolution equations then read 
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where K = N/12, and p = \{S — In R) is to be substi- 
tuted. The constraint equations are 
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where w stands for either Hort), and T ww denotes the 
semiclassical energy fluxes (whose explicit form is not 
needed here). We have set here /j = (a trivial solution 
of the scalar field equation fi. uv = 0), as we are deal- 
ing here with the evaporation of the BH rather than its 
formation. 

It is useful to re-express the system of evolution equa- 
tions @ in its standard form, in which R ruv and S ruv are 
explicitly given in terms of lower-order derivatives: 
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This form makes it obvious that the evolution equations 
become singular at R = K and R = 0. 



B. homogenous set-up 

In Ref. Q it was established that the spacetime inside 
a macroscopic evaporating CGHS black hole is (locally) 
approximately homogeneous. Motivated by this obser- 
vation, we shall restrict our attention to homogeneous 
solutions of the field equations. Namely, we consider so- 
lutions which only depend on t = v + u. In the homoge- 
neous framework the evolution equations (jj) become 



A. The model and field equations 

The semiclassical CGHS model [3] describes gravity in 
I + I dimensions coupled to a dilaton (f>. It also includes 
a cosmological constant A 2 and a large number TV >> I 
of identical scalar matter fields fi . We express the metric 
in double-null coordinates u, v, namely ds 2 = —e 2p dudv. 
The action is 

i / dudv[e~ 2 * (-2p, uv +4 )U 0,„ -\ 2 e 2 P) ~ \ (I) 
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where an over-dot denotes differentiation with respect to 
t. The constraint equations <j3j> now turn into a single 
equation 



The last term represents the semiclassical correction. For 
notational simplicity we set A = 1 [this choice amounts 
to a trivial shift p — >• p + ln(A)]. 



R - R3 + T = 0, 



C. The R = K Singularity 

The semiclassical evolution equations are singular at 
R = and R — K. In a black-hole interior R de- 
creases monotonically, starting from a macroscopic value 
(R>> K) at the horizon. Therefore, of these two singu- 
lar values, it is R = K which is first encountered. Since 
the evolution beyond this singularity is apriory unknown, 
it is not clear whether the semiclassical R = singularity 
will ever be encountered. For this reason we concentrate 
here on the R = K singularity (which is guaranteed to 
develop) . 

As it turns out, this singularity is characterized by the 

divergence of R and S, while R and S are finite. The 
evolution equations are dominated near the singularity 
by! 



R = — 



R 2 



2(R-K) 

The general solution is 

R (t) = K-\ 
S(*) = ~|t-«o| 
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(9) 



It depends on four arbitrary parameters B, B 2 , B^,to as 
required. We find that R and S are indeed continuous at 

t = to, but R and S diverge as \t — to| -1 / 3 . 

To simplify the description of the asymptotic behavior 
near R = K we re-define our variables as 



R=R 



Z = R + KS. 



The above near-singularity evolution equations then de- 
couple and take the simple form 
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D. 



Effective Hamiltonian for the near-singularity 
region 



Next we construct an effective Hamiltonian which 
yields the leading-order asymptotic behavior near R = 
K, with the aim of subsequently quantizing this dynam- 
ical system. 

For the i?-motion one finds the Hamiltonian [5f 



R 



(11) 



(up to an arbitrary multiplicative constant), where P is 
the momentum conjugate to R. For Z we obviously have 
the free-particle Hamiltonian Hz — |-P§, where P% is 
the momentum conjugate to Z. The overall leading-order 
effective Hamiltonian is H + Hz- 



Note that the effective Hamiltonian (fTTj) is independent 
of t, hence H must be conserved along any semiclassical 
solution R(t). Indeed one finds that R = dH/dp = 2P/R 
and therefore Eq. © yields a constant H, 
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III. 



QUANTIZATION OF THE SYSTEM AND 
EIGEN-FUNCTIONS 



We turn now to quantize our dynamical system, pro- 
ceeding as outlined in Sec. [I] Since the near-singularity 
dynamics of z is trivial, we shall focus here on the quan- 
tization of R (the trivial quantum evolution of z will be 
briefly addressed at the end of Sec. IV]) . The dynamics of 
R is thus described by a wave- function ^(R, t), which ad- 
mits the usual probabilistic interpretation of Schrodinger 
theory. The wave function evolves in time according to 
the Schrodinger equation 



ihd t ^ = HV. 



(13) 



Here H denotes the Hamiltonian operator, namely an 
operator version of the semiclassical Hamiltonian (Jill) , 
in which P is replaced by the momentum operator P = 
—ihdR. 

Since P is a differential operator, the ordering of P and 
the factor 1/ R matters, and different orderings may lead 
to different quantum theories. There are many possible 
orderings, and a priori it is not clear which is the "right" 
one. Here we choose the simplest symmetric ordering: 
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In 1+1 dimensions, the combination HG/c 3 turns out 
to be dimcnsionless (for any choice of length, time, and 
mass units). Since we have already chosen Gencral- 
Relativistic units c = G = 1, the constant h is dimcn- 
sionless. In what follows we choose the value h = 1. A 
different choice of ft will amount to a rescaling of the " En- 
ergy" parameter E (which we shortly introduce) — and 
correspondingly a rescaling of the time scale for quan- 
tum evolution. [12( 

The system's eigenstates i/je(R) are determined by 
the time-independent Schrodinger equation HtPe(R) = 
Eijje(R), namely 



did, 
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(15) 



The general solution of this equation is 

i, E (R) = C A Ai'(^E^ 3 R) + C B Bi'(-E 1 / 3 R) (16) 

(for E 7^ 0), where Ai and Bi respectively denote the 
Airy functions of the first and second kinds, the prime 



denotes a derivative with respect to their argument, and 
Ca,Cb are arbitrary coefficients. [13[ 

Note that both functions Ai(x) and Bi{x) are real, 
and are analytic at x = 0. It is remarkable that unlike 
the semiclassical singularity at R = — and despite 
the singular 1/R factor which explicitly appears in the 
Hamiltonian — all eigenfunctions i^e (R) are absolutely 
regular at R = 0. This fact is crucial to our approach: It 
implies that a generic wave packet ifj(R,t) propagating 
towards small R will not experience any singularity on 
approaching R = 0, hence it will admit a unique and 
well-defined evolution afterward. 

We may only use eigenfunctions which are square- 
integrable. Note the far-zone asymptotic behavior of the 
Airy functions: For x » 1 we have 
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and on the negative axis: 
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(17) 
(18) 

(19) 
(20) 



At the side of positive x (negative ER), the Bi eigen- 
functions diverge exponentially in |i?| 3 / 2 , and are hence 
non-normalizablc. On the other hand, the functions 
Ai {x) admit a well-behaved asymptotic behavior for 
both signs of R (exponential decay at one side and os- 
cillations at the other side). Thus, the permitted eigen- 
functions are 



i> E {R)=Ai{-E 1 / z R), 



(21) 



yielding an unbounded, non-degenerate, continuous spec- 
trum. 

In the next section we shall construct a well-localized 
initial wave packet far from the singularity [14J|, and the 
far-zone asymptotic behavior of the eigenfunctions will 
play a key role in the construction. We introduce new 
variables related to R and E, in order to give this asymp- 
totic behavior a more standard form of harmonic func- 
tions. 

First we define r(R) (and its inverse) to be 

r^^- sl gn{R)\R\V\ R=t^sign{r)\r\ 2 '\ (22) 

We also dchne a variable fc satisfying 

1 



E=-k 2 . 
2 



(23) 



We shall only use E > eigenfunctions for the con- 
struction of our wave packet (see next section). In terms 



of the new variables fc, r, the E > eigenfunctions 

ip k (r) = ip E{k) (R(r)) read 

ipk (r) = Ai [sign(r) | (3/2)fc r | 2/3 ] . (24) 

Their far-zone asymptotic behavior becomes: 

^(r)^- l(3/2) ^L r| cos(|fc|r + 7r/4) (\k\r » 1) 

v/tt 

(25) 

and 
V> fc (r)«-^^ — exp(-\kr\) (\k\r « -1) . 



2V^ 



(26) 



IV. 



CONSTRUCTING THE INITIAL WAVE 
PACKET 



Next we construct a well-localized initial wave packet 
^o(r) = ^(r, t = 0) from the system's eigenfunctions. 
In the analogous free-particle problem, it is often found 
convenient to consider gaussian wave packets. Inspired 
by the harmonic form of the asymptotic behavior (|25[) 
(along with the relation E ex k 2 ), which resembles a free 
particle, we shall seek an approximately-gaussian initial 
wave function in our problem too. This gaussian is pre- 
sumably centered at some initial location r$ > 0, with 
width Sro, and it propagates from large to small r with 
a certain averaged wave-number fco < 0. The desired 
fiducial gaussian is thus 



e lkor e 






(27) 



We design our initial wave function to be sharply 
peaked in both r-spacc and fc-spacc. Correspondingly, 
we demand Sr << r Q and Sk << |/c | , where 5k = 
1/8tq. These inequalities (which may always be achieved 
by choosing sufficiently large r$ and |/co|) also imply 
r |fc | » 1. 

To accord with the above requirements, we take our 
initial wave function to be 



*oW= / c(k)i> k (r)dk, 



where 



c(fc) = pe-' lro ~ k e - Sr °~ k2 /\ 



(28) 
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k = k — fco, and (3 is a certain normalization factor (de- 
termined below). [15| Note that c(fc) is just the fc-th 
Fourier coefficient associated with the fiducial gaussian 
([2"T]). Figure Q] below illustrates (for the specific case 
r = 4, fc = -40, Sr = 1/2) that this ^a(R) has the 
desired well-localized gaussian-like shape. 

We now analytically demonstrate that Eqs. (|28I29|) . 
combined with the above inequalities, indeed yield an 



approximate Gaussian. Since our wave packet is pre- 
sumably sharply peaked in both r-space and fc-space, the 
inequality ro|fco| >> 1 ensures that |fc|r » 1 is satis- 
fied throughout the phase-space region which dominates 
^>o (?") ■ We may therefore approximate the eigenfunctions 
tpk(r) by their far-zone harmonic form (|25[) . Also, in the 
factor cos(|fc|r + 7r/4) therein we may replace |fc| by — fc to 
facilitate analytic integration. [161 ] Furthermore, in the 
prefactor |(3/2)fcr| 1 / 6 we may approximate k by fco. We 
obtain the approximate asymptotic form 
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where a = — 7r ±/2 (-§-| ^o I ) X / ^ - Substituting this approxi- 
mate ipk{r) in Eq. (|28| and integrating over k, one ob- 
tains 
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where 70 
gaussian at r 



_ (-l)^yAF 



<5rn 



a/3 . This describes the desired 



ro plus an undesired one at r ~ — ro- One 
should recall, however, that the harmonic approximation 
([30| for ^(r) only holds for positive r (with |fc|r >> 1). 
For negative r the eigenfunctions are instead exponen- 
tially damped (by the huge factor exp(\ko\ro) >> 1), as 
seen in Eq. (|26|) . Therefore we are only left with the 
desired positive-r (approximate |l7[) gaussian: 



7o|r| 1/6 e <fc ° r 
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The normalization condition J | 'Q'o I («iR/<i7") <ir 

implies 70 = 
(i-H)V?n7 
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2 1 /12y^ w l/4: 
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1 now 
which in turn yields /3 = 



2 5/12, I .l/4| fco |l/6 • 

Summarizing this section, our initial wave function 
^o(r) is given by Eq. (J28J) combined with Eqs. (|29|) 
and ([M)l - which yield the desired well-localized, 
approximately-gaussian shape (|32[) . Switching back from 
r to R we obtain 

e ir S-^S a /4 i4i '[_J^i ii ] djfcj (33) 



*o(i?)=/3 



where, recall, fc = fc — fco. 

Figure [1] displays both the exact ^(-R) of Eq. (|3"3")l and 
its approximate form (|32[) (with r translated back to i?) , 
for the set of parameters ro = 4, fco = —40, and <5ro = 
1/2. The two graphs are indistinguishable, indicating the 
validity of the above approximations. 



V. TIME EVOLUTION OF THE WAVE PACKET 

The time evolution of ^> is dictated by the Schrodingcr 
equation (|13[) . which endows each component ipk in Eq. 
([28} with an extra phase factor exp(— iEt): 
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c(fc)Vfc(r)e-^ fc2)t dfc. (34) 
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FIG. 1: The initial wave function |\&o| as a function of R, 
constructed from Eq. (|33[) (solid black), and compared with 
the approximation ((32} (dashed gray). Note the non-existence 
of the "mirror" Gaussian on the negative side of R. Note also 
that ro = 4 corresponds to initial location R = 3 2 ' 3 2 « 4.16 . 



More explicitly, transforming r back to R, 

J— 00 " 

(35) 
We shall first inspect ^>(R,t) numerically, and then ex- 
plore it analytically. 



Numerical evaluation 

We numerically evaluated Eq. (|3"5} for the particular 
case ro = 4, fco = —40, 8rg = 1/2 (yielding 8k = l/8ro = 
2). 

To this end we approximate the integral over fc by 
a summation over a discrete set of equally-spaced val- 
ues of fc, separated by a small interval e, in the range 
fco — Afc>fc>fco + Afc. To achieve numerical pre- 
cision the parameters s, Afc must satisfy e « 8k and 
e -(Ak/2Sk) 2 <<x We choose here £ = (Sfc/20 = 0.1 and 
Afc = hSk = 10. The plots were insensitive to any further 
increase of Afc or decrease of e. [18l | 

Figure [2] shows a three-dimensional plot of ^(R,t). 
It clearly demonstrates how the incidenting gaussian 
bounces off R = 0. The time evolution roughly con- 
sists of three stages: (1) propagation from positive R 
toward R = 0, (2) the intermediate stage, near the mo- 
ment of bounce, and (3) propagation of the wave packet 
back from the neighborhood of R — towards large posi- 
tive R. On approaching the singularity (stage 2) the wave 
packet develops wiggles, but subsequently it bounces and 
becomes nicely-peaked again (stage 3). This behavior is 
seen more clearly in Fig. [3J which displays a few snap- 
shots of *(i?,i). 

We also verified numerically that the wave function 
preserves its unit norm upon evolution as expected. |19j 




FIG. 2: A three-dimensional plot describing the time evolu- 
tion of the wave function |$(it,t)|. The bounce at R = is 
evident. 



A. Analytic study 

Suppose that at some moment the wave function is 
located away from r = 0, in a range of r satisfying 
r |fco| >> 1- Then, as was discussed in the previous sec- 
tion, the oscillatory far-zone asymptotic behavior applies 
and we may represent the eigenfunctions by their ap- 
proximate form (|30p . It is useful to define an auxiliary 
function ^(r, t) by substituting this approximation for 
tpk(r) inEq. §Q: 



§(r, t) = a |r| 1/6 / c{k) e"^ 2 )* C os(-fcr + n/4)dk. 

J — oo 

(36) 
This function provides valuable insight into the dynam- 
ics of ^ as we shortly demonstrate. The integration is 
straightforward [in fact this integral is exactly the same 
as the time evolution of a free Schrodinger particle with a 
double-gaussian initial state as in Eq. (f3"Tj) ]. We present 
^ as 



■t*n- 



e"" u ■ e Sr <*r 



where ", 



and 
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r c (t) = r Q + k t, 



(37) 
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$_ and $ + arc certain phase factors (|$±| = 1) whose 
explicit form will not concern us. Evidently \l/(r, t) is a 
superposition of two gaussians which propagate at oppo- 
site directions, one located at r = r c (t) and the other at 
r = —r c (t). Both gaussians approach r = at time 
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FIG. 3: Plots of |\t| as a function of R at several differ- 
ent times, as it propagates towards the singularity (a), hits 
the singularity (b), and then recedes from it (c). The solid- 
gray lines and dashed-black lines represent the numerical 
calculations and analytic approximations, respectively. Fig- 
ure 3a displays three snapshots which, going from right to 
left, correspond to times t = 0, 0.04, 0.08. Fig. 3b Displays 
1*1 at t = thit = 0.1. Fig. 3c again displays three snap- 
shots which, going from left to right, correspond to times 
t = 0.12, 0.16, 0.2. The analytic approximation (dashed- 
black line) in the two right-most graphs of both figures 3a 
and 3c corresponds to Eq. (|41[) . In Fig. 3b the analytic ap- 
proximation corresponds to Eq. (|42[) . 
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(40) 



As long as the two gaussians comprising iff are remote 
from r = — namely, \r c (t)\ >> 5r(t) — the inequality 
|rfco| >> 1 is satisfied throughout both gaussians. Re- 
peating the argument which led to Eq. (|32l) above, we 
find that the true wave-function * will only contain the 



positivc-r Gaussian: 



— t^ — 



|*(r,t)| ~7|r| 1/b e <*■« 



|r c (t)| »<&■(*)• 

(41) 

This single-gaussian shape characterizes the wave func- 
tion in both stages 1 and 3 (namely, the propagation to- 
ward R = 0, and back from R = towards large positive 
R, respectively). It is illustrated in the two right- most 
graphs in both Figures [3^ and |3fc. These graphs dis- 
play both the exact wave-function (|35|) and the (visually- 
indistinguishable) corresponding single-gaussian approx- 
imation (|4ip . Notice that the gaussians in Fig. [3J: are 
wider than their counterpart in Fig. [3Ji, which simply 
expresses the monotonic growth of Sr(t). 

The intermediate stage 2 occurs at i ~ thiti when the 
gaussian arrives at the very neighborhood of r = 0. Dur- 
ing this stage |\&| develops numerous wiggles, as demon- 
strated in Fig. [3Jd (and also, to some extent, in the left- 
most graphs of Figs. [3b. and [3b). Qualitatively this be- 
havior may be interpreted as an interference pattern re- 
sulting from the overlap of the two gaussians comprising 
iff. In particular the right wings of both gaussians super- 
pose at a range of positive r values of order a few times 
dr. This superposition leads to an oscillatory interfer- 
ence pattern in "P, with a typical wave-number kg. Since 
|fco|5r >> 1, the far-zone harmonic approximation ([30)) 
holds at r ~ Sr > (for both gaussians), hence ^ rj ^ 
in this range. [2(| 

This qualitative picture only yields a rough approxi- 
mation for ^ in stage 2, because near r = the far- 
zone approximation breaks. However, one can show (by 
matching two different asymptotic approximations) that 
at the very moment of bounce the following approxima- 
tion holds: 

*(r,t = t hu ) = Ae-^Vh(fc 2 /2) 1 / 3 J R(r)]. (42) 

where q = (Sr^ — 2ito) and A is a certain coefficient. 
Figure [3h> displays both this approximate expression and 
the full function ^(r, t = thit)i demonstrating a nice 
agreement. 

Finally, we comment on the quantum dynamics of the 
other degree of freedom Z . Since Z evolves as a free par- 
ticle [Sec Eq. (|10[) ] the quantum formulation adds noth- 
ing to its semiclassical dynamics: The wave-function ty z, 
which describes the quantum dynamics of Z ', will just 
propagate freely, totally ignoring the R = singularity. 



VI. SEMICLASSICAL EVOLUTION BEYOND 
THE SINGULARITY 

In both stages 1 and 3, the wave- function is sharply 
peaked at 



r = \r c (t)\ = k \t-thit\- 



(43) 



back from r to R, we obtain 

R(t) = B \t - t \^ , 



(44) 



with Bq = |3/co/v / 2| 2 ^ 3 and to = thiti m mu consistency 
with Eq. ©. 

One of our main goals was to obtain the semiclassical 
extension beyond the R = singularity. This extension is 
now determined (within the near-i? = approximation) 
by the last formula: It is just a time-reflection of the 
semiclassical solution at t < to- 

Note that the parameters Bq , to which characterize the 
extension at t > to are fully determined by the semiclas- 
sical evolution at t < to. In particular, in Ref. [5J it was 
found that for a macroscopic (M >> K) CGHS black 
hole Bo ~ [(3/2)V^M] 2 / 3 , where M denotes the (re- 
maining) black-hole mass. 

The other degree of freedom Z just evolves steadily at 
t ~ to (at the leading order) ; Namely, 



Z(t) ~ z +zi(t-t ), 



(45) 



where Zq, Z\ are constants which are again determined by 
the semiclassical evolution at t < to. 

Note that the extension (|4"4"|) . which satisfies B(t > 
to) = B(t < to), is the only one which conserves H across 
the singularity, as may be seen from Eq. (fT2|) (which holds 
at both t < t and t > to). 

It should be emphasized that the exact time-reflection 
symmetry exhibited in Eq. ([44| only characterizes the 
leading-order semiclassical evolution near R = 0. Beyond 
the leading order R and Z do couple [as may be seen, for 
example, by substituting S = (Z — R)/K in Eq. flSJ], 
which spoils the time- reflection symmetry of R(t). The 
semiclassical evolution at t > to beyond the leading order 
may be obtained by taking Eq. (j4"4")) and the above linear 
expression for Z(t), translated back to R and S, as initial 
conditions for the full system (|5|6p . 



By the correspondence principle, this orbit should corre- 
spond to some semiclassical solution. Indeed, translating 



VII. DISCUSSION 

It is a common belief that Quantum Gravity will even- 
tually cure the spacetime singularities, which are known 
to occur in classical and semiclassical gravitational col- 
lapse (as well as in cosmology). Here we demonstrated 
this idea by applying a simplified quantization proce- 
dure to the spacelike R = singularity inside a two- 
dimensional CGHS [3j evaporating BH. This allows us 
to explore the extension of spacetime physics beyond the 
singularity. 

Our quantization procedure is based on the minisu- 
perspacc approach, in which all degrees of freedom as- 
sociated with inhomogeneities are erased. The system's 
evolution is then described by a wave-function, whose 
time evolution is dictated by a Schrodinger-like equation, 
which wc solve analytically. 

Our most crucial observation is that quantum evo- 
lution is perfectly regular (in fact, analytic) across the 



semiclassical R = singularity. (This is far from obvi- 
ous, because the Hamiltonian operator itself is singular 
at R = 0.) This analyticity gives rise to unique quan- 
tum evolution across the singularity. We find that an 
initially-localized wave packet spreads and fluctuates on 
approaching R = 0, but then it bounces back toward 
positive R values. This is demonstrated in Figs. [2] and [3] 
Following the bounce, the wave packet becomes sharply- 
peaked again, and approaches a well-defined semiclassical 
orbit. 

The main objective of this paper was to explore the 
extension of spacetime physics beyond the R = sin- 
gularity inside the BH. The wave-packet dynamics de- 
scribed above indicates that the strong quantum fluctu- 
ations which develop at the very neighborhood of R = 
quickly die out, giving rise to a new semiclassical phase 
beyond the singularity. The quantum dynamics also de- 
termines the structure of spacetime (namely, the space- 
time dependence of the two semiclassical variables R and 
Z) in this new patch, by dictating the boundary con- 
ditions for R and Z across the singularity: First, both 
variables are continuous. Second, R bounces there at 
a (locally) time-symmetric manner, whereas Z evolves 
steadily across the singularity. These boundary data for 
both R, Z and their time derivative (which in fact amount 
to specifying the parameters B^, to, Zq, z\ of Sec. VI) give 
rise to a unique and well-defined semiclassical Cauchy 
evolution at the future side of the spacelike singularity, 
the region denoted by D in Fig. |4] 




FIG. 4: Spacetime diagram of a CGHS black hole which 
forms by gravitational collapse and subsequently evaporates. 
The thick null line at the bottom left represents the collapsing 
massive shell. The event horizon is marked by a dotted line. D 
represents the future Cauchy development of the singularity 
line R — 0. The domain marked by J is discussed in the text. 

Interestingly, the situation found here resembles the re- 
cent developments in Loop Quantum Cosmology (LQC) 
in many respects. Several analyses showed that in var- 
ious homogeneous cosmological models (open or closed 



universe, with or without cosmological constant) the clas- 
sical big-bangsingularity is replaced in LQG by a regu- 
lar bounce. [9( The same behavior (quantum singular- 
ity resolution, and a consequent bounce) was obtained 
here in a two-dimensional BH — though within a more 
simplified quantum formulation. There appears to be 
an important difference, though: In our case the very 
epoch of bounce is marked by strong quantum fluctua- 
tions, whereas the cosmological bounce in LQC is typi- 
cally free of such strong fluctuations. 



An obvious extension of this research is to obtain the 
specific structure of semiclassical spacetime [namely, the 
functions R(u, v) and S(u, v)] beyond the singularity. In 
the present work we obtained the required boundary data 
for these functions at the future side of the singularity. 
The calculation of post-singular R(u, v) and S(u, v) may 
naturally be divided into three stages: (i) The asymp- 
totic behavior in the very neighborhood of the singular- 
ity is already given in Eqs. (|4"4")l and P5|) . It applies as 
long as R << K. [21| (ii) In an intermediate stage of 
larger t — to values, these simple asymptotic expressions 
no longer hold. Furthermore, R and Z no longer evolve 
independently. Still (as long as t — to < < M/K) the ho- 
mogeneous approximation applies, and the evolution of 
R(t) and S(t) is determined by the ordinary differential 
equations (|5I6[) . (iii) For even larger values of t — to, the 
large-scale inhomogencous character of spacetime comes 
into expression, and the evolution of R(u, v) and S(u, v) 
must be determined by evolving the partial differential 
equations ((4]). In both stages (ii) and (iii), the evolution 
can be analyzed either numerically or by certain ana- 
lytic approximations (which are beyond the scope of the 
present paper) . A numerical simulation of the semiclas- 
sical CGHS field equations is underway [IflLLlI], an d per- 
haps it will be possible to extend this numerics to the 
post-singularity region as well. By this way it should be 
possible to determine the semiclassical spacetime [partic- 
ularly the functions R(u, v) and S(u, v)] in the entire fu- 
ture "domain of dependence" of the spacelike line R = 0, 
the triangle-like region denoted by D [which is actually 
D + (R = 0)} in Fig. HI 



Note, however, that after the structure of the domain 
D + {R = 0) is exposed, one will still face the most chal- 
lenging stage in the task of composing the overall post- 
singularity semiclassical spacetime: Revealing spacetime 
physics in the "wedge" between D + (R = 0) and the ex- 
ternal asymptotically-flat region (the domain denoted J 
in Fig. 0|). A priory it is not even clear if this "piece" 
is connected (as suggested by the ATV point of view |8| , 
which assumes a fiducial Minkowskian !ft 2 manifold) or 
disconnected (as may be suggested by the more common 
paradigm, see e.g. Fig. 5 in Ref. pj). It seems that new 
concepts and insights will be needed in order to explore 
this remaining piece of the puzzle. 
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